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■ ABSTRACT 

> 

lO . We discuss morphological properties of the large-scale jets in powerful radio 

sources, which are now observed at radio, optical and X-ray frequencies, in order 



O , to determine the origin of their X-ray radiation and the nature of the particle 

Q ! acceleration responsible for their multiwavelength emission. We show that mod- 

^ ! eling knots of the these objects as stationary regions of energy dissipation within 

^! the uniform and continuous jet flow leads to several problems. Such problems are 

2 ! especially pronounced if the observed X-ray emission is due to inverse-Compton 

'. scattering of the CMB radiation. We explore another possibility, namely that 

^ ' the knots represent moving and separate portions of the jet matter, with excess 

kinetic power. We suggest a possible connection between this scenario and the 
^ \ idea of intermittent /highly modulated jet activity. The proposed model can ex- 

5^ \ plain some morphological properties of quasar jets — like high knot-to-interknot 

brightness contrasts, frequency-independent knot profiles and almost universal 
extents of the knot regions — independently of the exact emission mechanism 
responsible for producing the X-rays. In this context, we consider different pos- 
sibilities for the production of such X-ray radiation, and discuss the related issue 
of particle acceleration. We conclude that the appropriate process cannot be 
designated for certain yet. We suggest, however, that X-ray observations of the 
jet in 3C 120 (in which intermittent jet activity is quite distinct) already seem 
to support the synchrotron origin of the X-ray emission, at least in this object. 



Subject headings: acceleration of particles — radiation mechanisms: nonthermal — 
galaxies: jets — X-rays: general 
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1. Introduction 

Growing numbers of large-scale jets in powerful extragalactic radio sources are now 
being observed at optical and X-ray frequencies (see, e.g., Stawarz 2004). New, and in many 
ways unexpected, results call for a wide discussion of the physical conditions within these 
objects, and in particular of the processes responsible for accelerating radiating particles 
to the required ultrarelativistic energies. Until now, few observations had gone beyond the 
integrated spectral properties of multiwavelength emission from large-scale jets. However, the 
high angular resolution of Hubble Space Telescope (HST) and Chandra X-ray Observatory 
(CXO) have given us the opportunity to analyze the morphological characteristics of jets at 
optical and X-ray photon energies, allowing for the comparison of jet morphology at different 
frequencies. Here we discuss the structure of large-scale quasar jets, in order to determine 
the origin of their X-ray radiation and the nature of the particle acceleration responsible for 
their multiwavelength emission. 

CXO has detected aboTit ten quasar jets with observed X-ray luminosities Lx ~ 10^^ — 
]^q45 gj,g g-i g^^^ projected lengths ranging from a few tens up to a few hundred kpc^. The 
X-ray emission is usually much too intense to be explained in terms of synchrotron radiation 
by the same, single-power-law electron energy distribution that produces the radio-to-optical 
emission. Nor can it be explained in terms of synchrotron-self Compton (SSC) radiation of a 
homogeneous, few-kpc-long knot region, unless the energy density is far above the equiparti- 
tion value, which would imply extremely high jet powers (e.g., Chartas et al. 2000; Schwartz 
et al. 2000, for the case of PKS 0637-752). This led to the idea that the X-ray emission results 
from inverse-Compton up-scattering of Cosmic Microwave Background (CMB) photons — 
'external inverse-Compton' (ETC) emission — by the low-energy part of the electron energy 
distribution (Tavecchio et al. 2000; Celotti et al. 2001; Harris & Krawczynski 2002). Such a 
model, however, requires highly relativistic jet velocities at distances of a few hundred kpc 
from the active nuclei (bulk Lorentz factors F ~ 10), in apparent disagreement with radio 
observations (e.g., Wardle & Aaron 1997). On the other hand, the ETC model was claimed 
to be additionally supported by energetic considerations (GhiseUini & Celotti 2001). In this 
respect, discoveries of large-scale X-ray jets in the Gigahertz-peaked (GPS) objects PKS 
1127-145 and B2 0738+313 (Siemiginowska et al. 2002, 2003a, respectively), as well as in the 
high-redshift source GB 1508+5714 {z = 4.3, Siemiginowska et al. 2003b) are particularly 
interesting. X-ray emission of large-scale quasar jets, if due to the EIC process, has also 
cosmological implications, as pointed out by Schwartz (2002). Therefore, it is important to 
discuss whether Comptonization of the CMB within highly relativistic flow is really the most 



see the web page http://hea-www.harvard.edu/XJETS by D. Harris; see also Stawarz (2004) 
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plausible explanation of the CXO quasar jet observations, and to look for possibilities of dis- 
criminating between this and other models (e.g., Aharonian 2002; Dermer & Atoyan 2002; 
Stawarz & Ostrowski 2002). Analysis of the morphology of the emitting regions constitutes 
an interesting approach to this problem. 

The most apparent characteristic of quasar jets is their knotty morphology with high 
knot-to-int erknot brightness contrast, but also with distinct (in some cases) intcr-knot diffuse 
emission. In addition, knot profiles seem to be frequency- independent, and knot extents are 
similar when observed at radio, optical and X-ray photon energies. Detailed observations of 
the 3C 273 jet (Jester et al. 2001) reveal also that spectral changes along the fiow arc not 
correlated with brightness changes — this may be a general characteristic of these objects. 
Finally, in some cases spatial offsets between the maxima of the X-ray and radio emission 
within the knot regions were noted (e.g., PKS 1127-145, Siemiginowska et al. 2002). It is 
not clear whether all of these features can be explained in a framework of models involving 
extended shock waves within continuous jet flow. In fact, we argue that the morphological 
characteristics cannot be explained in this way, and that substantial modifications of the 
standard picture are required. Such modifications are especially needed if the X-ray emission 
of quasar jets is due to the EIC process. Wc propose that at least some aspects of the 
HST and CXO observations can be understood in terms of intermittent (modulated) jet 
activity. In this context, we comment on both the external- Compton scenario of quasar jet 
X-ray emission, and on a model of boundary layer acceleration and resulting high-energy 
radiation (Stawarz & Ostrowski 2002). In particular, in § 2 we emphasize the problems with 
modeling X-ray knots of quasar jets as stationary regions of particle acceleration. In § 3 we 
consider the possibility that the knots are moving sources of non-thermal radiation, propose 
a possible connection of this scenario to models of intermittent jet activity, and discuss the 
particle acceleration processes possibly involved in such a scenario. The discussion and final 
conclusions are presented in §§ 4 and 5, respectively. 



2. Stationary knots 

In the most common version of the EIC model, in which the knots are identified with 
strong stationary shocks through which the jet matter fiows continuously, none of the mor- 
phological features mentioned in the introduction are straightforwardly expected. The knots 
appear to be localized along the jet, yet stationary shocks are likely to be highly oblique 
(with respect to the jet axis) and significantly elongated in the flow direction, if the jet is 
confined by external pressure (Sanders 1983; Komissarov & Falle 1997). Large-scale jets in 
powerful radio sources display approximately cylindrical geometry, and therefore it is thought 
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that they must be confined. The nature of their confinement is unknown, but it has been 
shown that the low-density and high-pressure nonthermal medium of the radio lobe can be 
an efficient confinement agent (Begelman & Cioffi 1989). If we are not seeing the shape of 
the shocks themselves, could the emissivity patterns give rise to the knot morphology? It 
seems unlikely, since the cooling time of low-energy electrons responsible for up-scattering 
CMB photons to the observed X-ray energies is extremely long when compared with the 
lifetime of electrons emitting synchrotron optical photons in an equipartition magnetic field. 
In addition, the EIC emission is not very sensitive to the magnetic field structure or slow 
adiabatic losses within the uniform, well-coUimated fiow. Therefore, in the case of the EIC 
model and a continuous jet with knot regions marking the positions of the stationary (e.g., 
reconfinement) shocks, one should rather expect extended X-ray knot emission, with ap- 
proximately constant flux along the jet (if the jet bulk velocity is roughly constant on the 
large scales), and with knot proflles reflecting directly the large differences in radiative and 
adiabatic cooling times of electrons radiating in radio, optical and X-rays. This is in dis- 
agreement with observations, and hence modifications to the EIC model are required. One 
possibility is to assume an inhomogeneous structure of the emission regions, consisting of 
a number of small, far-from-equipartition and adiabatically expanding clumps of radiating 
plasma (Tavecchio et al. 2003). 



2.1. EIC clumping model 

The EIC model was claimed to be strongly supported by energetic arguments (Ghisellini 
& Celotti 2001). Let us briefly re-discuss this issue in the context of the clumping scenario. 
If the jet is composed of a magnetic field of unknown intensity B (measured in the comoving 
frame^), protons, and radiating electrons, one can find a value of the magnetic field that 
minimizes the total kinetic power of the emission regions for a given observed synchrotron 
luminosity Lgyn, observed volume of the emitting plasma V, ratio of proton-to-electron energy 
densities 1] = u'^/u'^, spectral energy distribution of the radiating ultrarelativistic electrons, 
bulk Lorentz factor F, and jet inclination 9. This is because the total jet power is Ltot = 
Lb + Le + Lp, where the Poynting flux is Lb oc B^, the bulk kinetic power carried by the 
radiating electrons is Le oc u'^, and the bulk kinetic power due to the protons is Lp oc Wp. 
As the electron energy density for a given total synchrotron luminosity is u'^ oc B~^, one has 



^Hereafter we use primes to denote quantities measured in the jet comoving frame, while unprimed 
quantities refer to the observer rest frame. However, we do not prime the magnetic field induction B or 
electron Lorentz factor 7, noting instead that they always refer to the emitting plasma rest frame. 
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B=B, 



= 0, 



(1) 



one can find that the magnetic field Bcr minimizing Ltot is the one for which Lb = Lc + 
Lp. In other words, one gets condition for the cquipartition between Ponting fiux and the 
jet kinetic power carried by the particles (see Appendix A). It was shown (Ghisellini & 
Celotti 2001) that this minimum power condition is consistent with X-ray observations of 
the large-scale jet associated with the quasar PKS 0637-752, if its X-ray flux is due to the 
EIC emission of a homogeneous knot region with radius R ~ 10^^ cm and F ~ 10. One 
should emphasize, however, that this agreement is obtained by invoking certain assumptions 
about the value of 77 and details of the electron energy distribution. In addition, allowing 
for inhomogeneous structure of the emission region changes the energetic conclusions. In 
particular, the magnetic field that minimizes the total power for the uniform knot model, 
Bci, docs not coincide with the minimum energy field evaluated in the framework of the 
clumping scenario for the radiating plasma, Bcr- Instead, one has 



where / is the filling factor of the clumps (see Appendix A). Hence, if the magnetic field 
estimated from the observed ratio of inverse- Compton to synchrotron luminosities agrees with 
the critical magnetic field computed for the homogeneous knot model, B ~ B^r, then any 
strong clumping of the jet (i.e. / -C 1) will cause large departures from the minimum power 
condition for the radiating plasma, B <^ Bcr- In addition, in order to discuss energetics of the 
whole clumpy jet, one has to make some assumptions about the non-radiating matter between 
the clumps, which cannot be verified by observations. As a result, energetic arguments 
supporting the EIC model cannot be applied to the clumping scenario. 

Losing a reliable diagnostic of jet energetics is the price one has to pay for explaining a 
rapid decrease of the X-ray fiux outside the knots. In the expanding, radiating plasma, the 
magnetic field is expected to decrease in proportion to some power of the emitting region 
radius, B oc r~"^, depending on the assumed magnetic field structure. Therefore, one has 



where Ldc and Lgyn are, respectively, the observed inverse- Compton and synchrotron lumi- 
nosities, u'^^^ is the energy density of the CMB radiation in the comoving frame of the jet 
(which is constant along the flow if the bulk velocity of the expanding region is constant), 
and u'q oc B^ is the energy density of the jet magnetic held. Hence, knot proflles in the 



Bcr^r'^'B, 



cr 



(2) 
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EIC model with adiabatic expansion should always be different when observed at X-ray and 
radio frequencies (both produced by the low-energy part of the electron spectrum), unless 
very rapid expansion of the radiating plasma, i.e., strong clumping of the jet, is involved. In 
this case, however, the contribution from synchrotron self-Compton (SSC) emission cannot 
be ignored. Note that the main argument in favor of the EIC emission model over the SSC 
model was that the latter requires clumps and departures from equipartition, in order to 
explain the strong X-ray emission (Schwartz et al. 2000). The second problem attributed to 
the ElC-clumping model was that rapid expansion of the emitting clumps would inevitably 
create a large population of very low energy electrons. They should in turn produce op- 
tical emission via the EIC process, possibly exceeding the synchrotron radiation (and the 
observational limits) at these frequencies. 

Energy losses of the low-energy electrons, radiating via synchrotron radiation in the 
radio band and via EIC emission in the X-ray domain, are dominated by adiabatic cooling, 
d'y/dr — —A^jr. This process leads to the decrease of the electron energy with linear size 
of the emitting region r according to 

7 (X (4) 

and also to the decrease of the inverse-Compton luminosity at given observed frequency 
according to 

LeieOcr-^(^-^), (5) 

where s is the spectral index of the electron energy distribution (Appendix B). For s = 2.5 
and two-dimensional expansion (i.e., A = 2/3), one gets a linear decrease of the observed 
X-ray luminosity with r. In such a case, in order to reduce the X-ray flux by, say, one 
order of magnitude, one has to allow all the emitting clumps to expand within the single 
knot region by, at least, one order of magnitude in their linear sizes. Thus, electrons up- 
scattering CMB photons to the observed < kcV energies, 7 < (^^x/^cmb)^^^ ~ 1000 5"^, 
will end up at energies 7 < 200 after such expansion, giving rise to the observed EIC 
emission at e < 50 eV (see Tavecchio et al. 2003). This possibility, although not completely 
excluded, contradicts the synchrotron origin of optical emission from quasar jets, which has 
been confirmed repeatedly by polarization measurements. Hence, the EIC model requires 
a very sharp cut-off in the injected electron energy distribution and/or relatively limited 
expansion of the clumps — both of which call for fine-tuning of the model parameters — in 
order not to overproduce optical inverse-Compton radiation. 

Small sizes of the emitting plasma in the clumpy jet can also amplify the energy density 
of the synchrotron photons to a level leading to domination of the SSC component over the 
EIC one. Note that, for the electron energy spectrum usually considered in the EIC model, 
the photon energy range of the SSC emission corresponds roughly with the spectral range of 
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the EIC radiation. The ratio of observed SSC to EIC luminosity is 

-^ssc 1 f S \ ^ f Lsyn \ f Rc 

Leie ^ N{l + zY \10J ^1044erg/sy' VlOpc 

where A'" is the number of the emitting clumps with radius Rc per knot, z is a redshift of the 
source and S — l/r(l — is the jet Doppler factor for its bulk velocity /3 = (1 — F"^)^/^ and 
inclination 9 = cos~^ /i (Appendix C). It is always possible to adjust the unknown parameters 
S, Rc and N in order to obtain Lssc > L^ic- As a result, it is possible to model the observed 

X-ray emission of large-scale quasar jets by SSC emission, and hence there are no obvious 
constraints on the large-scale jet parameters from X-ray observations. In particular, the 
Doppler factor of the jet in the SSC model can be smaller than that required by the EIC 
model. 




2.2. Synchrotron X-rays? 

Modifications to the EIC model, which are required in the case of knots treated as 
stationary regions of energy dissipation, decrease the model's attractiveness. Therefore, 
alternative scenarios postulating a synchrotron origin of the X-ray emission should be con- 
sidered. However, in this case the X-ray photons would have to be produced by a different 
electron population than the one responsible for the radio-to-optical emission, as the X-ray 
flux is usually above the extrapolated radio-to-optical power-law continuum; or else devia- 
tions from a single power-law behavior in the electron energy distribution — high energy 
particle pile-up in particular — have to be invoked. Furthermore, to reconcile the observed 
flat X-ray spectral indices with strong cooling of the highly energetic electrons and large 
observed sizes of the emission regions (sec discussion in Aharonian 2002). one must have 
particle acceleration occurring over extended volumes. This can be illustrated by rewrit- 
ing the observed propagation length of the electrons radiating synchrotron photons with 
frequency i/gyn in a form 



-'"^'^^^ i + io-3Pi^_i(i + .)4 Kw^J ^P^' 

where t[^^ is the radiative loss time scale, B^/^ = i?/10~^ G, invcrsc-Compton radiative losses 
on the CMB are included and we put 5 ~ F as appropriate for small jet viewing angles 6 
(Appendix D)^. 



^Let us mention that, as emphasized by Jester et al. (2001), the apparent discrepancy between this spatial 
scale for the optical synchrotron emission and the observed extent of the optical jet in 3C 273 cannot be 
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The turbulent second-order Fermi acceleration process provides a plausible way of ob- 
taining high-energy particle pile-up, if the time scale for particle escape from the acceleration 
region is much longer than the radiative loss time scale (other possibilities will be discussed 
in the later sections). This is because such a process favors the particle energy for which 
the acceleration time scale is equal to the radiative loss time scale. With steady injection of 
low-energy particles, turbulent acceleration piles up radiating particles at this critical energy 
As discussed by Ostrowski (2000) and Stawarz & Ostrowski (2002), boundary shear layers of 
large-scale jets are ideal sites for such a process, as the inevitable flow instabihties very likely 
provide the required amount of turbulence, and a sheared magnetic field configuration does 
not enable radiating electrons to escape easily from the acceleration region. For conventional 
assumptions about the jet energetics (including energy cquipartition between the magnetic 
field and radiating electrons), the time scales estimated by Stawarz & Ostrowski (2002) lead 
to sufficient pile-up of the high-energy electrons to produce synchrotron X-rays with a lu- 
minosity comparable to the observed values, without fine-tuning of the jet parameters. In 
addition, the continuous nature of such an acceleration process can explain the absence of 
cooling effects in the observed optical and X-ray continuum. 

Boundary layer acceleration, as initially proposed to explain some aspects of the X-ray 
jet emission observed by CXO, cannot solely account for the morphological characteristics 
of quasar jets, and their knotty nature in particular. In this model the lifetime of the X-ray 
emitting electrons is extremely short, and hence the emitting region sizes are related to the 
extent of the particle acceleration sites and not to the spatial scales of the electron coofing, 
as in the EIC and SSC scenarios. Also, the interplay between localized shock acceleration 
and continuous, turbulent energization of the radiating electrons can result in offsets between 
the radio and X-ray emitting regions, while in the EIC and SSC models the maximum of the 
radio emission should have the same location within the knot as the maximum of the X-ray 
radiation. As mentioned in the introduction, spatial offsets between radio and X-ray knot 
maxima are observed in some large-scale quasar jets. However, a simple model of turbulent 
acceleration does not explain directly why the jets — being continuous flows — exhibit high 
knot-to-inter knot brightness contrasts, although one can imagine some possible solutions. 
For example, it was suggested (Stawarz & Ostrowski 2002) that the knots can form by the 
shock compression of the energetic electrons previously accelerated at the boundary layer 
upstream of the shock, in accordance with the discussion by Begelman & Kirk (1990). 



resolved solely by invoking a sub-equipartition magnetic field and/or highly relativistic velocities. 
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2.3. Spectral constraints 



Detailed spectral studies at optical frequencies could allow us to distinguish between an 
inverse-Compton and a synchrotron origin for the X-ray emission from large-scale quasar jets. 
For example, the presence of spectral cut-offs at infrared-to-optical photon energies would 
argue strongly for the EIC model, while spectral flattenings would support synchrotron 
models involving high-energy particle pile-ups. However, optical emission detected from a 
few quasar jets^, which is most probably synchrotron in nature as indicated by polarization 
measurements, reveal a variety of spectral shapes, and hence no definitive statements can 
be made yet. In particular, optical spectral indices of the knots in powerful jets range from 
Oo ~ 0.5 in the brightest part of the 3C 273 jet (Jester et al. 2001) up to > 1.6 in 3C 279 
(Cheung 2002). Interestingly, spectral hardening at ultraviolet frequencies was noted for the 
3C 273 jet (Jester et al. 2002). These observations indicate the presence of 7 ~ 10^ — 10^ 
electrons producing synchrotron infrared-to-optical photons, for an assumed equipartition 
magnetic field B ~ 10~^ G, but do not give tight constraints on their spectral distribution. 
The broad-band energy spectrum of the radiating electrons is usually modeled in the EIC 
scenario as a single power- law, n'^ij) oc 7"*, with spectral index s ~ 2.5 and sharp cut-offs 
at 7min ~ 10^ and 7inax ~ 10^ — 10^ (e.g., Sambruna et al. 2002). Except for its simplicity, 
however, such an assumption on a form of ^^(7) has neither theoretical nor observational 
justification. Moreover, it cannot be excluded that the electron energy distribution n'^i'j) 
extends up to very high energies, 7 ~ 10^, as seems to be the case in a few kpc-scale FR I 
jets (see discussion in Stawarz 2004). Another point is that the broad-band radio-to- X-ray 
spectra for some of the knots in the large-scale jets associated with quasars 3C 273 and PKS 
1136-135 cannot be explained by the EIC model, indicating with little doubt the synchrotron 
origin of the X-rays (Marshall et al. 2001; Sambruna et al. 2002). Therefore, existing optical 
spectra of the knots in quasar jets cannot yet exclude different scenarios for modeling the 
radio-to- X-ray continua of these objects. 



3. Moving knots 

As we have discussed above, modeling knots of large-scale quasar jets as stationary 
regions of energy dissipation within the uniform and continuous jet flow leads to several 
problems in explaining the observed X-ray emission, especially in the framework of the EIC 
scenario. Below we explore another possibility, namely that the knots represent portions 
of a relativisticaly moving jet with excess kinetic power. Such a situation is expected if 



see the web page http://home.fnal.gov/~jester/optjets by S. Jester; see also Stawarz (2004) 
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the activity of central engines in radio-loud AGNs is intermittent, or highly modulated. 
The frequency-independence of the knot sizes, as well as high knot-to-interknot brightness 
contrasts, can be understood if the radiating particles are present mainly within the separate 
portions of the jet flow — the knots — as compared with the interknot regions. This is true 
almost independently of the emission process involved, and in particular it applies also to 
the EIC model. In addition, in moving knot scenarios the brightness changes along the jet 
path do not have to correlate with the spectral changes, if they are controlled by distinct 
physical processes, i.e., if the former are due to changes in the kinetic energy of the flow 
resulting from intermittent or modulated jet activity, while the latter are determined by 
particle acceleration and coohng processes. These facts make the idea of intermittent jet 
activity promising for explaining the morphologies of CXO quasar jets. 



3.1. Modulated jet activity 

The idea of intermittent jet activity was presented in the context of 'partial' radio jets 
associated with the powerful radio galaxies 3C 219 (Bridle et al. 1986; Clarke et al. 1992) 
and 3C 288 (Bridle et al. 1989). It was also proposed that the restarting jet model can, in 
a natural way, reconcile the presence of the extended radio emission observed in a few GPS 
sources (e.g., 0108-1-388) with the inferred young ages of GPS compact radio structures, < 10^ 
yrs (Baum et al. 1990; Stanghellini et al. 1990). We note that extended radio emission around 
GPS objects can sometimes reach Mpc scales, as in the cases of B1144 -1-352 (Schoenmakers 
et al. 1999) and B1245-I-676 (Marecki et al. 2003a). In addition, as discussed by Reynolds & 
Begelman (1997), intermittent jet activity can explain the general overabundance of compact 
radio objects when compared with the number of classical (extended) doubles. In particular, 
Reynolds & Begelman (1997) showed that the observed distribution of sizes for extragalactic 
radio sources can be explained in the framework of a simple evolution scenario if one assumes 
that the jet-like active period persists only over ~ 10^ years, and refreshes on the time scale 
of > 10^ years. Interestingly, similar time scales can be obtained in the thermal-viscous 
instability model for the accretion disk (Siemiginowska & Elvis 1997; Janiuk ct al. 2003). 
This suggests that the jet activity — and hence the whole evolution of the radio source — 
could be driven by unstable accretion disks surrounding supermassive black holes (see in this 
context Marecki et al. 2003b). 

Non-steady jet activity, albeit on longer time scales than discussed above, was also 
suggested in order to explain some morphological features observed in the so-called Giant 
Radio Galaxies (linear sizes of the order of Mpc; Subrahmanyan et al. 1996). In a few sources 
belonging to this class, intermittent jet activity was confirmed directly by observations of 
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double-double radio structure (Schoenmakers et al. 2000a). In particular, observations of 
B1834+620 (Schoenmakers et al. 2000b), J0116-473 (Saripalli et al. 2002), and B1545-321 
(Saripalli et al. 2003) show a pair of 'inner' radio jets which propagate within a diffuse radio- 
emitting medium formed by previous jet activity, and which display FR II radio morphology. 
Such inner radio structures are characterized by lower radio luminosity, larger axis ratios of 
the lobes, and weaker hot-spots than classical doubles with similar sizes, in agreement with 
theoretical predictions for restarting bipolar jets (Clarke & Burns 1991; Kaiser et al. 2000). 

A restarting jet propagates not within the unperturbed intergalactic or intracluster 
medium of the active galaxy, but within the interior of the radio lobe created by the previous 
jet activity (e.g., Begelnian & Cioffi 1989). The number of double-double radio sources is 
small, probably because this kind of morphology can be formed only if the new pair of jets 
is created not earher than 10^ — 10^ yrs after the old ones, in order to allow for an increase 
in the density of the cold gas in the relict radio cocoon (Kaiser et al. 2000). Otherwise, the 
rarefied and mostly nonthermal medium formed during the previous period of jet activity 
prevents formation of new lobes and hot-spots^. For the quasar jets considered here, we 
require the time scale between the two subsequent jet active epochs to be not longer than 
10^ yrs. This is because the linear sizes of the knot regions are consistent with a 10"^-year 
duration of jet activity, while distances between subsequent knots are on average consistent 
with quiescent periods an order of magnitude longer. Therefore, we characterize the quasar 
jets in terms of modulation of the jet kinetic power, rather than in terms of restarting flows. 
This idea is strongly supported by radio observations of such structures (see, e.g.. Bridle et 
al. 1994). The respective time scales for activity and quiescence can differ from object to 
object, but on average they seem to be in agreement with the ones postulated above. In 
some sources, however, they can be smaller than average, resulting in a more continuous 
flow (as in 3C 273?). 

We also note that Galactic X-ray transients are known to produce radio emission which 
is widely modeled in terms of multiple jet-like ejection events (see review by Fender 2003). 
Modulated jet activity of the large-scale jets considered here constitutes an interesting link 
between quasar and microquasar phenomena, although the physical processes controlling 
intermittent /highly modulated activity of the central engines in both types of sources can 
be different. As a result, the time scales for the active and quiescent periods in galactic and 
extragalactic jet sources do not have to scale in a simple way with the mass of the central 
compact object. 



^On the other hand, the time scale of the quiescent periods in double-double sources cannot be longer 
than 10^ yrs, because after this time the old lobe structure is not expected to remain visible (Komissarov & 
Gubanov 1994). 



3.2. Energy dissipation and peirticle acceleration 



The moving portions of the jet have to be sites of particle acceleration. Tavecchio et al. 
(2003) rejected this possibility, arguing that either continuous or episodic particle accelera- 
tion acting within moving knots cannot reproduce the required electron energy distribution. 
However, as the synchrotron spectrum at frequencies higher than the radio band is unknown 
in most large-scale quasar jets, such a statement is premature (see § 2.3). Note that Tavec- 
chio et al. (2003) base their conclusions on the assumed electron energy distribution. In 
reality, many different processes can control the spectral behavior of radiating electrons, and 
it is not certain — from both observational and theoretical points of view — how well a 
single power-law describes the spectrum. In the framework of the intermittent jet activity 
model considered in this section, turbulent particle acceleration processes due to interaction 
of the jet with the surrounding medium are expected to lead to the formation of a piled-up 
distribution of energetic electrons radiating synchrotron X-rays. The formation of internal 
shocks within the relativistic, highly modulated jet flow is also very probable. In the latter 
case, one can expect the EIC process to dominate production of the X-rays, whereas the de- 
ominant radiation mechanism is less clear in the former. Realistically, the situation can be 
quite complex, involving an ensemble of multiple shocks (formed due to both velocity irreg- 
ularities and propagation instabilities of the flow) as well as the interplay between coherent 
and stochastic particle acceleration processes. 



3.2.1. Internal shocks 

The formation of shocks in extragalactic jets is most likely related to intrinsic veloc- 
ity irregularities of the flow, caused by the non-steady operation of the central engine, as 
proposed by Rees (1978). Such internal shocks are not stationary, but instead are moving 
together with the radiating matter, possibly with highly relativistic bulk velocities. Internal 
shocks are typically characterized by a low efficiency of converting the beam kinetic energy 
to internal energy of radiating particles (see, e.g., Lazzati et al. 1999; Spada et al. 2001). 
This is a very convenient situation for the objects considered here, because their observed ra- 
diative luminosities, if corrected for Doppler beaming effects, are several orders of magnitude 
smaller than the estimated jet kinetic powers (Lj ~ 10^^ erg s~^), and hence the fraction of 
the kinetic energy flux of the jet dissipated to the radiating particles can be small. 

Let us consider the simplest situation, in which two colliding portions of the jet flow 
differ only in bulk velocity. Then, for respective bulk Lorentz factors r2 > Fi ^ 1, a double 
shock structure (symmetric in the comoving frame of the contact discontinuity for the same 
densities pi = P2) develops. The bulk Lorentz factor of the contact discontinuity, which is 
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also the Lorentz factor of the shocked radiating plasma, is then 

Tct ~ Vt^2. (8) 

Knowing solely Fi and r2 (and assuming an equation of state for the downstream matter), 
one can estimate velocities of the forward and reverse shocks from standard shock-jump 
conditions (Blandford & McKee 1976, Appendix E). The extent of the shocked region along 
the jet, in the contact comoving frame denoted by primes, is 

Al' ^2c\(3'JAt', (9) 

where \f3'g-^\ is the velocity of the forward and reverse shock fronts and At' is the time since 
formation of the double-shock structure. Observations of the knots at the maximum syn- 
chrotron frequency z/max can give some constraints on the maximum value for A/'. This is 
because the time interval At' cannot be larger than the radiative loss time scale for the most 
energetic electrons radiating at t'max, in order to ensure activity of the whole shocked region 
at this frequency. Therefore, for a small viewing angle of the jet, the observed extent of the 
shocked region can be limited according to 

A/ < Ux = 2 C Tet Cd('^max), (10) 

where t'^.^^ = t'.ad(^max) is the radiative loss time scale for a given synchrotron frequency 
^max (Appendix D). Figures 1-3 show the parameter /max as a function of the magnetic field 
strength in the radiating plasma, for z/max = 10^^ Hz and different bulk Lorentz factors Fi 
and F2. Note that, for a wide range of bulk Lorentz factors of the colliding portions of the 
jet, the extent of the double-shock structure in the observer frame cannot be larger than 
~ 1 kpc, independent of the magnetic field intensity. This value, interestingly, is consistent 
with the observed lengths of the optical knots in quasar jets. Therefore, if dynamical time 
scales controlling the development of the double-shock structure ensure that A/ < /max for 
the maximum synchrotron frequency occurring in the optical band, then the ElC/internal 
shock scenario can explain frequency-independent knot profiles. 

However, there are two important problems with the model discussed above. The first 
is that radio polarization of the quasar jets indicates a magnetic field which is on average 
parallel to the jet axis. As an example we mention the jet in PKS 0637-752, where the X-ray 
emission as observed by CXO is strong only in the regions characterized by a longitudinal 
magnetic field (Schwartz et al. 2000). Such a situation is inconsistent with the internal shock 
model, because compression of the magnetic field should result in a transverse orientation 
(Laing 1980). The second problem is the fine-tuning required between the dynamical time 
scales controlling development of the double-shock structure, and the radiative loss time 
scales of the radiating particles, in order to obtain almost universal extents of the optical/X- 
ray knots, ~ 1 kpc. 
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3.2.2. Multiple oblique shocks and turbulent acceleration 

We now discuss particle acceleration processes related to ensembles of shocks which 
are not comoving with the radiating matter, i.e., have a different pattern velocity than 
the jet. Numerical simulations reveal the presence of such a structure at the jet boundary 
(e.g., Norman et al. 1982), although it is possible for these shocks to encompass the whole 
jet width. The particle acceleration process related to the ensemble of many shocks can 
be different from single shock acceleration. This issue was considered in the hterature in 
different contexts starting from Blandford & Ostriker (1980), and discussed further by Spruit 
(1988); Achterberg (1990); Schneider (1993); Melrose & Pope (1993); Pope & Melrose (1994); 
Melrose & Crouch (1997); Marcowith & Kirk (1999) and Gieseler & Jones (2000). Physically, 
what distinguishes between single and multiple shock acceleration is that in the latter case, 
particles can change their momenta between subsequent shock events due to adiabatic and 
radiative losses, and can also escape from the regions between of the shocks, thus avoiding 
further acceleration. As shown by Schneider (1993) and Marcowith & Kirk (1999), the 
interplay between time scales representing radiative cooling, particle escape and advection 
can result in the formation of different particle energy spectra, with spectral flattenings, 
curvatures, cut-offs or high-energy bumps (spectral pile-ups). 

Let us consider the appropriate time scales for multiple shock acceleration in the con- 
text of the knots in quasar jets. Following Schneider (1993), we assume that the time scale 
for acceleration and subsequent decompression at each shock is much shorter than the time 
scales connected with radiative losses, t'^^^^, and particle escape from the region between sub- 
sequent shocks (i.e., the region where the nonthermal radiation is produced), t'^^^, as well as 
the interval between the subsequent shock acceleration events, t'^. Then, the average parti- 
cle energy distribution resulting from multiple (roughly identical) shock acceleration events 
depends on the ratios of these three characteristic time scales. Note first, that if particle 
escape is controlled by diffusion from the region between the shocks, which is characterized 
by a thickness D, then the condition t'^^ > 4ad reads 

L>>0.04 5Zf^ pc (11) 

for strong turbulence, where the particle mean free path is comparable to the particle gy- 
roradius, and for inverse- Compton losses comparable to synchrotron losses (Appendix F). 
Thus, if D is of order of the jet radius (~ 1 kpc), then particle escape is inefficient compared 
to the radiative losses, independent of particle energy. Also, if the mean separation of the 
shocks is less than ~ 1 kpc, particles cannot escape before passing through many shocks, 
and therefore they experience multiple (and not single) shock acceleration. In such a case, 
as shown by the analysis of Schneider (1993), one should expect the formation of a fiat 
power- law spectrum for particles satisfying t'^.^^ > t[j, a steep power-law spectrum (but not 
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necessarily an exponential cut-off) for particles satisfying t[g^ < t\^, and a spectral pile-up at 
particle energies for which t'^.^^ ~ t[j. 

Now consider what the shock separation has to be in order to obtain a piled-up distri- 
bution of particles emitting the observed synchrotron keV photons. If the radiating plasma 
between the shocks is highly relativistic (with the respective bulk Lorentz factor F » 1), 
then the condition ^ad ~ reads 



where d' is the shock separation as measured in the radiating plasma rest frame (Appendices 
D and F). The parameter d' as a function of F, for the observed critical synchrotron photon 
energy hucr ~ 1 keV and different magnetic field intensities, is shown on figure 4. For 
F ~ 3—10 and B ~ 10^^ — 10"'^ G one gets d' ~ 10—100 pc. For larger shock separations, pile- 
ups arc expected to occur at lower energies of the synchrotron photons. However, these values 
should be considered as illustrative only, because the analysis above refers to a very simplified 
situation, in which we neglect the effects of turbulent second-order Fermi acceleration and 
assume identical properties for all the shocks. For reasonable large-scale jet parameters, 
however, an ensemble of shock waves present within the extended knot region can result 
in the formation of a flat power-law electron energy distribution followed by a high-energy 
pile-up announced by its X-ray synchrotron emission. This situation is similar to the model 
discussed by Stawarz & Ostrowski (2002), but here the time scale between subsequent shock 
events, and not the second-order Fermi acceleration time scale, tpjj, determines the critical 
electron energy. The ratio of these two time scales, for B_4 ~ 1 and electron Lorentz factor 
7 ~ 10^, can be rewritten as 



where V\ = Pac is the Alfven speed in the radiating plasma (Appendix F). Therefore, we 
conclude that an ensemble of oblique shock waves within the knot regions of large-scale 
quasar jets can act in the analogous way to second-order Fermi acceleration, as discussed 
by Stawarz & Ostrowski (2002). Both processes, which in a realistic situation can be hard 
to separate, are due to the interaction of the propagating jet with the surrounding medium. 
Both of them can result (for reasonable jet parameters) in a pile-up of high-energy electrons 
radiating X-ray synchrotron photons, and both are consistent with the observed longitudinal 
magnetic field. 



d' ~ ct' 



(12) 




(13) 
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4. Discussion 

Inverse-Compton scattering of the CMB radiation field within a highly relativistic fiow 
is commonly believed to explain relatively strong X-ray emission detected by CXO from 
large-scale quasar jets. This model, as initially proposed, is attractive because of its sim- 
plicity and energetic efficiency. On the other hand, it is inconsistent with morphological 
properties of these objects, namely with the frequency-independent knot profiles, if the 
jet flow is continuous and the knots are stationary (e.g., if they are manifestations of the 
reconflnement shocks through which the jet matter is flowing at a constant rate). It was 
proposed (Tavecchio et al. 2003) that the inconsistency can be resolved by assuming that 
the radiating jet region consists of small, adiabaticaly expanding clumps. However, as dis- 
cussed in § 2.1, by invoking an inhomogeneous structure of the emitting regions, one loses 
the advantages of the EIC model over other models proposed in the literature. In particular, 
when clumps are incorporated into the EIC model, emitting regions generally have to be 
far from the minimum power condition. In this situation, it becomes much more difficult to 
arrange for the radiating particles to be confined to the required volumes. Also, knot profiles 
will be observed as almost frequency-independent only if the radiating regions expand very 
rapidly, starting from very compact clumps. Such a situation leads to the problem that 
contributions from the SSC process cannot be excluded anymore. Furthermore, fine-tuning 
of the parameters is required in order not to overproduce optical emission by the EIC process 
at the final stage of the clump evolution. One should emphasize that, if the X-rays from 
large-scale quasar jets — modeled as continuous flows with localized and stationary particle 
acceleration sites — are indeed inverse-Compton in origin, then clumping of the knot regions 
is unavoidable. In the case of the EIC scenario this is because strong adiabatic losses are 
required in order to obtain knots instead of uniform and continuous emission along the jet, 
while in the case of the SSC process it is because the energy density of synchrotron photons 
within the emitting regions is required to be larger than in the homogeneous case in order to 
provide the observed level of X-ray emission. Therefore, we consider the clumpy jet model 
(as proposed by Tavecchio et al. (2003)) as unlikely. 

We note that complex substructure of the knots is suggested by the time variability of 
the kpc-scale jet in the PR I radio galaxy M 87 (e.g., Harris et al. 1997, 2003). In PR I 
sources, the X-ray emission observed by CXO is most probably synchrotron in origin. Also, 
detailed studies of the Centaurus A jet radio and X-ray radiation seem to indicate a spine- 
boundary layer morphology (Hardcastle et al. 2003). This can be taken as an indication 
that dissipation processes acting at the jet boundary can be crucial for understanding the 
high-energy radiative output of extragalactic large-scale jets, and their X-ray emission in 
particular (Stawarz & Ostrowski 2002, § 2.2). We stress that the nature of the quasar jet 
clumps — as localized small sites of particle acceleration — is unknown. It was proposed 
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(Tavecchio et al. 2003) that the clumpy structure can form when the jet crosses a region 
filled with small clouds of ambient medium, or that the clumps can be related to magnetic 
reconnection events. 

Regardless of the specific mechanism of X-ray production, the morphological properties 
of quasar jets (inferred from multiwavelength observations) can be understood if the knots 
are not stationary regions of energy dissipation within the continuous jet flow, but instead 
are moving, separate portions of the jet matter. In § 3.1 we outlined the possible connection 
between this idea and highly modulated jet activity. In this case, a variety of particle accel- 
eration processes (including reconnection, acting in a way analogoTis to second-order Fermi 
acceleration) can play a role. Unfortunately, optical observations, which could be crucial 
in discriminating among models, cannot give definitive answers at the moment. Therefore, 
detailed spectral and morphological studies are still needed. Let us mention in this context 
one especially interesting object, which could be a very promising target for such multiwave- 
length studies: the large-scale jet associated with the Seyfert galaxy 3C 120. In this object, 
intermittent jet activity seems to take place, while the X-ray emission of the knot seems to 
be exceptionally inconsistent with the EIC model. 

3C 120 {z — 0.033) displays large-scale radio structure, with a one-sided bending jet 
(observed up to roughly 100 kpc from the core) and diffuse, amorphous radio emission 
extending in several directions up to > 400 kpc (Walker et al. 1987). VLBI observations 
reveal complex jet structure with superluminal motions of some jet components at the 1 — 100 
pc scale (e.g.. Walker ct al. 2001). This, together with the large bending of the flow at kpc- 
scales, indicate a small jet inclination to the line of sight. At ~ 2 — 3 kpc from the host 
galaxy, an elongated radio knot, oriented perpendicular to the flow axis, is observed. Long- 
term VLA monitoring (Walker 1997) indicates that this knot is stationary (or moving only 
subluminally) . Farther away the jet bends and fades, but is manifested again at ~ 20" from 
the core as a bright, extended radio knot, oriented obliquely with respect to the jet path. 
The radio spectral index ~ 0.65 remains constant along the flow, although brightness 
variations are large. Optical observations reveal polarized emission aligned with some parts 
of the radio jet in the flrst few kpc from the core, but with morphology not corresponding 
directly to the radio one (Hjorth et al. 1995). Surprisingly, the 20" knot turn out to possess 
an X-ray counterpart with a 0.2 — 2 keV luminosity ~ 10^^ erg s~^, as shown by ROSAT 
observations (Harris et al. 1999). Optical upper hmits (placed below the interpolated radio- 
to- X-ray continuum) allow one to estimate the optical-to- X-ray power-law slope for this 
region as cKq-x < 0.35. 

The EIC model applied to the X-ray knot of the 3C 120 jet gives an uncomfortably 
small viewing angle of the jet, 9 ~ 1.5°, and a very high bulk Lorentz factor F ~ 40 (Harris 
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& Krawczynski 2002). Note that, in such a case, the deprojected distance of the X-ray knot 
from the host galaxy is on the order of Mpc. On the other hand, radio observations together 
with optical upper limits exclude the possibility that the radio-to- X-ray continuum results 
from the synchrotron emission of a single power-law electron energy distribution. Also, 
if the X-ray emission is synchrotron radiation, for example due to an additional electron 
population as suggested by Harris et al. (1999), the fiat spectral index and large extent of 
the knot call for continuous particle energization^. We suggest that the X-ray emission of 
the knot in the 3C 120 jet is indeed synchrotron in origin, and is a manifestation of the 
pile-up of high-energy electrons by some particle acceleration process acting continuously 
within the radiating region. This region can be further identified with a separate portion 
of the jet matter, possibly moving with relativistic velocities (but not as high as required 
by the EIC model), ejected from the active nucleus some 10^ yrs ago. This is consistent 
with the morphology of the 3C 120 extended radio lobes, which also indicate intermittent 
jet activity. 



5. Conclusions 

We propose that the knots observed in large-scale quasar jets in radio, optical and X- 
rays represent portions of the jet with an excess of kinetic power due to intermittent high 
activity of the jet engine. The appropriate time scales for jet activity and quiescent epochs 
(~ 10^ yr and ~ 10^ yr, respectively) are roughly consistent with the ones estimated in 
different contexts for powerful radio sources (Reynolds & Begelman 1997; Siemiginowska & 
Elvis 1997). The proposed scenario can explain some morphological properties of quasar jets, 
independently of the exact emission process responsible for production of X-rays. Indeed, 
the dominant emission process cannot be designated for certain yet. It is possible that the 
X-ray emission of quasar jets is due to inverse-Compton scattering of the CMB by low- 
energy electrons accelerated in internal shocks. On the other hand, it is also possible that 
the observed keV photons are synchrotron in origin, being a manifestation of high-energy 
electrons piled-up at the critical energy ~ 100 TcV due to the second-order Fermi process 
and/or multiple shock acceleration. The fact that an exceptionally high Doppler factor is 
required to explain the 3C 120 X-ray data in the framework of the EIC scenario, suggests that 
synchrotron emission is likely to be the dominant process in this object, at least. Whatever 
emission mechanism prevails in most large-scale quasar jets, intermittent jet activity assures, 
in a natural way, high knot-to-interknot brightness contrasts, frequency-independent knot 



^Aharonian (2002) proposed an alternative explanation, involving synchrotron emission of very high 
energy protons. 
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profiles, and almost universal extents of the knots (~ 1 kpc). 
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A. Minimum power condition 

The total power of a cylindrical relativistic jet with bulk Lorentz factor F > 1 and 
radius R, consisting of magnetic field with comoving energy density u'q = B^/Stt, radiating 
ultrarelativistic electrons with the comoving energy density u'^ and protons with u'^ = r]u'^, 
is 

i-tot = Lb + Le + Lp (Al) 

where the Poynting fiux is 

Lb = TTi?^ cF^ u'ji (A2) 
and the power carried by electrons and protons is, respectively 

Le = ttR^ cF^ tig and Lp = r] L^. (A3) 

For a given electron energy distribution (i.e., the respective moments (7) and (7^)), the 
electron energy density can be related to synchrotron luminosity of some radiating volume 
V'by 



Hence 



2m,c2 it) B^r'iU 



For a radiating region modeled either as a moving or steady source, one has the transforma- 
tion to the observer frame 

T' 17" 

syn _ ^ ^syn /.^v 

v "53 V ^^^^ 

where the Doppler factor of the jet emission is 5 = 1/F(1 — Pn) and 9 = cos^^/x is the jet 
viewing angle (see, e.g., Sikora et al. 1997). Note that, for small jet viewing angles (5 ~ F), 
the Poynting fiux and the bulk power of radiating electrons do not behave in opposite ways 
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with respect to F and B, because Lb oc B"^ while Le oc F ^ B ^ (cf. Ghisellini & Celotti 
2001). Setting 

/f)T.. . \ 

= (A7) 



one can find that the magnetic field i?cr minimizing Ltot is the one for which Lb = + Lp. 
This value is given by 

Let us next consider a clumpy jet, consisting of identical clumps of radiating matter 
with comoving volumes VJ, which are present within the region of volume V. The clump 
filling factor can be defined as 

NV 

f - -yr- m 

Assuming that the jet velocity is uniform within the jet region under consideration, the 
Poynting flux for the radiating matter is then 

Lb ^ nR^ cT^ f u'b. (AlO) 

This follows from the fact that, given the magnetic energy of each clump ^ = V^, 
the total magnetic energy of the radiating matter flowing through the jet surface vri?^ is 
N ^/V = f u'b- Similarly, the bulk kinetic powers of the emitting electrons and protons 
(for an assumed 77 the same as in the homogeneous knot model) are, respectively, 

Le^rfcF^/iie and L^ ^ rj L^. (All) 

However, one has also 

K^^-j^ (A12) 

where L'^^^ ^ = L'^^^/N is the intrinsic synchrotron luminosity of each clump. As a result, 
the bulk kinetic energy of the particles computed for the radiating regions of the clumpy jet 
is the same as that computed for the homogeneous knot model. Hence, the magnetic fleld 
minimizing the total power of the radiating clumpy matter is 

4r = r'/'5cr. (A13) 



B. Adiabatic losses in the EIC model 

Adiabatic losses of ultrarelativistic electrons with Lorentz factor 7, within an expanding 
region of linear size r, obey 

^ = -^2 (Bl) 

dr r 
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(where A — 1 [2/3] for three-dimensional [two-dimensional] expansion), and lead to the 
evolution of the electron distribution according to 

A^'(7,r)= (^)^iV^(7o). (B2) 

Here A^'(7, r) = n'^{'y)V' , n'^{'y) oc 7~* is the electron energy distribution, Nq{'Jo) is the 
number of electrons injected impulsively at r = tq with initial energies 70, and 

T=(^)"\o (B3) 

(see, e.g., Atoyan & Aharonian 1999). 

The inverse- Compton luminosity due to electrons with given energy 7eic is 

[i/L,(r)]eic (X n;(7eic) V = N'ij,,,, r). (B4) 
Therefore, the adiabatic decrease of the ETC luminosity at a given frequency i^eic oc 7^;^ is 

[^L.(r)]eic _ iV'(7eic, r) _f r V^^'~'^ _ (35) 



[z/L,(ro)]eic A^o( /cic 



C. SSC emission in the clumping model 

The ratio of the observed SSC luminosity to the observed EIC luminosity is 

(e.g., Stawarz et al. 2003). The comoving energy density of the CMB radiation at redshift z 

4 

where oT^ ~ 4 x 10~^^ erg cm~^. The comoving energy density of the synchrotron emission 
in the clumpy jet, for clumps with radius i?c, is 

The intrinsic synchrotron luminosity of each clump, I/syn.c related to the observed syn- 
chrotron luminosity of the whole knot region, Lgyn, by 
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where N is a, number of identical clumps. Hence, one finally obtains 



Leic 47rcaT4 \1 + J 6^ N Rl' 



syn 



(C5) 



D. Radiative energy losses 

The time scale of radiative losses via synchrotron and inverse-Compton processes, where 
the latter is dominated by Comptonization of the CMB, as measured in the comoving frame 
of the radiating plasma is 

^ 1^ u'^{l + X) ^^'^ 

where X = i^cmb/'^B (see Appendix C). The electron Lorentz factor 7 is related to the 
observed frequency of the synchrotron emission, 1/, by 

3e 5B ^ , , 

i/ = 0.29 -7^ p2 

47rmeC (1 + ^) ^ ^ 

and therefore, for a given observed synchrotron frequency i/gyn, one obtains 

W ~ 1 + 3|Zr r2 5-2 (1 + ''syn • [^^> 

Due to inverse-Compton losses, for a given z/gyn the propagation length of the synchrotron- 
emitting electrons as measured in the observer rest frame, Zrad ~ ^^Kad^ cannot be arbitrary 
large, and passes through a maximum as a function of F or S (Jester et al. 2001). 



E. Internal shocks 

If the colliding portions of the jet flow differ only in bulk velocities {Pi and (32, re- 
spectively), one can find the bulk velocity of the contact discontinuity from the energy and 
momentum conservation laws: 

= A£i±|i£i (El) 

-L 1 + -L 2 

(e.g., Komissarov & Falle 1997). Hence, for highly relativistic velocities one has the contact 
discontinuity bulk Lorentz factor 



\/r\r^. (E2) 
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After the collision, a double-shock structure develops. This structure is symmetric in the 
comoving frame of the contact discontinuity (denoted by primes), and therefore the extent 
of the shocked region along the jet is 

Al' ^2c\f3'JAt' (E3) 

where is the shock velocity in the contact frame and At' is the time since the collision. 
In order to find the shock velocity, one has to solve the shock-jump conditions between 
the upstream plasma, moving with velocity and the downstream plasma, moving with 
velocity f3ct- This can be done conveniently in the comoving frame of the upstream plasma 
(denoted by (1)), to obtain the Lorentz factor of the shock 



r(i) 

^ sh 



(rL'^ + 1) [7ct (r2^ - 1) + 1]2 ^^^^ 



\ 7ct (2 - 7ct) (rS^ - 1) + 2 



(Blandford & McKee 1976). The constant 7ct can be identified with the ratio of specific 
heats of the upstream matter (we assume 7ct = 5/3). Next, the shock velocity in the 
contact frame can be found by transforming to the downstream plasma frame. For large 
Fi and F2, is expected to be subrelativistic. 

The shocked region is a moving radiating portion of the jet flow, and therefore its 
deprojected observed length is 

Al = 5ct Al' (E5) 

With the requirement At' < t'^^, where is the lifetime of the radiatively cooling electrons 
in the downstream plasma moving with bulk Lorentz factor Fct (Appendix D), one can write 

^<2c|/?:j4d (E6) 

Oct 

to obtain for small jet inclinations (i.e., for 5ct ~ F^) 

AZ<2c|/3^jFet4d- (E7) 



F. Multiple shocks 

The time scale for particles to escape from the region of the shock ensemble can be 
estimated as 

4c ~ — , (Fl) 
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where D is the thickness of the considered region and k ~ Ae(7) c/3 is the diffusion coefficient 
of the electrons with mean free path Ae(7) and Lorentz factor 7. Assuming that Ae(7) is 
comparable to the electron gyroradius, 

Ae(7)-r, = ^, (F2) 



one obtains 



^^1^^ (1 + X) (F3) 

where the radiative loss time scale 4ad given in Appendix D. 

The time scale between subsequent shocks can be estimated as 

t'^-^ (F4) 

where d' is the mean separation of the shocks in the comoving frame of the radiating plasma 
between the shocks (moving with bulk Lorentz factor F » 1). Second-order Fermi accelera- 
tion can be characterized by the time scale 

where the velocity of the scattering centers involved in the turbulent acceleration Vgc can 
be identified with the Alfven speed, Va, and the particle gyroradius can be taken for Ae(7) 
(Stawarz & Ostrowski 2002). Hence, one obtains the ratio 

^--(3l-^,d'B(3lr' (F6) 

''FII UleC 

where Pa = Va/c. 
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Fig. 1. — Maximum extent of the emission region in the internal shock model for the syn- 
chrotron frequency t'max = 10^^ Hz, as a function of the magnetic field B (at redshift z = 0.6). 
The solid line corresponds to the parameters Fi = 3, F2 = 9, Fct = 5.16 and /^g^ = 0.1826. 
The dashed line corresponds to the parameters Fi = 5, F2 = 15, Fct = 8.64 and P'^^i = 0.18. 
The dotted hne corresponds to the parameters Fi = 7, F2 = 21, F^ = 12.11 and (3'^^ = 0.1793. 
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Fig. 2. — Maximum extent of the emission region in the internal shock model for the syn- 
chrotron frequency t'max = 10^^ Hz, as a function of the magnetic field B (at redshift z = 0.6). 
The solid line corresponds to the parameters Fi = 3, F2 = 6, Fct = 4.23 and (3!,^ = 0.1179. 
The dashed line corresponds to the parameters Fi = 5, F2 = 10, Fct = 7.06 and P!,^ = 0.116. 
The dotted hne corresponds to the parameters Fi = 7, F2 = 14, F^ = 9.89 and (3'^^ ~ 0.115. 
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Fig. 3. — Maximum extent of the emission region in the internal shock model for the syn- 
chrotron frequency i/max = 10^^ Hz, as a function of the magnetic field B (at redshift z = 0.6). 
The sohd hne corresponds to the parameters Fi = 3, F2 = 3.3, Fct = 3.15 and P'^^ ~ 0.0167. 
The dashed hne corresponds to the parameters Fi = 5, F2 = 5.5, Fct = 5.24 and = 0.0162. 
The dotted line corresponds to the parameters Fi = 7, F2 = 7.7, F^ = 7.34 and (31^ — 0.016. 
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Fig. 4. — Separation of the shocks in the muhiple shock model corresponding to the critical 
synchrotron photon energy hu^- = 1 keV, as a function of bulk Lorentz factor F (for a redshift 
z = 0.6). The solid line corresponds to the magnetic field B = 10~^ G, the dashed line to 
B = 5- 10-5 G and the dotted line to B = 10"^ G. 



